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The position-representation wave function for multi-photon states and its equation of motion are 
introduced. A major strength of the theory is that it describes the complete evolution (including 
polarization and entanglement) of multi-photon states propagating through inhomogeneous media. 
As a demonstration of the two-photon wave function's use, we show how two photons in an orbital- 
angular-momentum entangled state decohere upon propagation through a turbulent atmosphere. 
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There are two approaches to solving problems in quan- 
tum mechanics: quantum field theory (QFT), and wave 
mechanics (WM). In WM, the fundamental physical enti- 
ties are "particles" , whose collective state is described by 
a wave function. To treat few-particle systems in WM, 
such as the helium atom with two electrons, one formu- 
lates and solves the two-electron wave equation in posi- 
tion space. In QFT, the fundamental physical entities are 
fields, which are decomposable into modes, each of which 
can have various numbers of excitations. It is clear that 
there is a one-to-one correspondence between modes in 
QFT, and states in WM. 

Each approach has its realm of preferred applicabil- 
ity. For example, one does not usually treat the helium 
atom with QFT, but one does use this theory when treat- 
ing high-energy electron collision experiments. Just as 
one most often uses WM in atomic physics, we advo- 
cate the use of this approach for few-photon phenom- 
ena, such as those encountered in elementary quantum 
information schemes (quantum cryptography, one-way 
quantum computing, and linear optics quantum compu- 
tation). However, as far as we know, a complete photon- 
wave-mechanics (PWM) theory of electromagnetism has 
not been introduced. Development of the PWM descrip- 
tion of multiple photons leads toward completion of the 
quantum description of electromagnetism, which must in- 
clude entanglement and dccoherence not yet treated. 

In this paper, we briefly review the one-photon wave 
function in coordinate space, of which there are several 
proposed 0, 0, 1, i, S, i, 0. Choosing one of these, 
we introduce a two-photon wave function, along with its 
tensor equation of motion, which we call the two-photon 
Maxwell-Dirac equation. Its relationships to other well- 
known formulations, such as quantum electrodynamics 
(QED) of few-photon wave packets [8|, and vector-field 
classical coherence theory Q are then brought to light. 
These connections suggest that we have chosen the most 
appropriate single-photon formalism upon which to base 
a generalization to multiple photons. Choosing a dif- 
ferent single-photon formalism does not give these close 
relations and leads to different equations of motion, nor- 



malization conditions, and Hilbert-space scalar products. 
We also demonstrate the use of the two-photon wave 
function in a calculation of quantum-state disentangle- 
ment for a pair of spatially-entangled photons traveling 
through the atmosphere. The two-photon formalism is 
readily extended to multiple-photon states. The the- 
ory gives a lucid view of the photon as a particle-like 
quantum object, making a pedagogical link between stan- 
dard quantum wave mechanics and quantum field theory, 
which must be equivalent according to current under- 
standing. 

Much of the confusion surrounding the definition of 
the photon wave function in the position representation 
arises from the non-localizability of the photon and the 
corresponding absence of a position operator [13, [H|. 
This is due to the fact that the photon has zero mass, and 
is a spin-1 object, with only two independent spin degrees 
of freedom. Of the coordinate-space single-photon wave 
functions proposed [E 0j S 0> S B > we nn d the most 
useful is the Birula-Sipe formulation,defined in terms of 
the localization of photon energy @, H, 0] 
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rather than, for example, the Landau-Peierls non-local 
number density wave function [l[. Here <jffl is a three 
component, complex vector labeled by a — ±1 for posi- 
tive(negative) helicity 0,0). In vacuum this wave func- 
tion satisfies the Dirac-like equation [1, H, 0) [H> 



ihdt$ {1) = H <F (1) = ftcE 3 V x 



(2) 



and the zero-divergence condition V • = 0. Here E3 
is a Pauli-like matrix that changes the sign of the neg- 
ative helicity component in Eq.(l). The curl and diver- 
gence operators are understood to act on the upper and 
lower components of separately, h is Planck's con- 
stant (which cancels in Eq.([2])), and c is the speed of light 
in vacuum. The single-photon Hamiltonian is ftcSsVx. 
Equation ^ and the zero-divergence condition are for- 
mally equivalent to the classical Maxwell equations, as 
can be seen by substituting 
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where D^ + ) and B( + ) are the positive-frequency parts of 
the electric-displacement and magnetic-induction fields, 
and eo (//o) is the vacuum permittivity (permeablity). 

In a linear, isotropic medium, eo and po in Eq.Q are 
replaced by spatially dependent functions e(x) and p(x), 
and the equation of motion Eq.([2]), is changed by the 
material interaction to [2j 

ihdt^ = H = (V + VL) x *W , (4) 

where the speed of light is constructed from the local 
values of permittivity and permeability in the medium, 
u(x) = l/y / e(x)/x(x), and the matrix L has the form, 

L(x) = (lln Ve(x)/i(x) + S x In Ve(x)//i(x)) /2. (5) 

Here 1 is the identity matrix and is a Pauli-like ma- 
trix that interchanges the two helicity components of 
Eq. ([1} . In a medium the photon Hamiltonian is given by 
hvT,z (V + VL) x. The modified divergence condition in 
a medium is (V + V-L) ■ \?W = 0. 

The integrated square modulus of the photon wave 
function over all space gives the expectation value of the 
photon's energy 



*W (x,t) t *W (x,i)d 3 :r = (£7i) 



(6) 



One can associate with this wave function a local proba- 
bility density p(x,t) = (x,t) j (x,i) / (Ei), and 
current density j(x,i) = (x, i) 1 * sfW (x, t) / (Ei), 
that obey the continuity equation [2], Here s is a 
vector composed of the three spin-1 matrices. This prob- 
ability density and current density are denned in relation 
to the photon energy, not photon number. 

The appropriate scalar product is best formulated in 
momentum space, where a local photon-number proba- 
bility density is well defined [ijj]. Transformed into the 
position representation, the scalar product is found to be 
a nonlocal integral, consistent with the absence of a local 
photon particle-density amplitude [3] . The fact that pho- 
ton wave functions representing orthogonal states are not 
orthogonal with respect to an integral of the form Eq. © 
is consistent with the well-known nonexistence of local- 
ized, orthogonal spatio-temporal modes in QED Q. 

We now propose that the two-photon wave function 
'J'*- 2 ) (xi,X2,t), which is related to the probability am- 
plitude for finding the energies of two photons localized 
at two different spatial positions xi and X2, at the same 
time t, with the photons in any polarization state, can 
be constructed from single-photon wave functions as 

VP ( Xl , X 2 , t)=Y, Clm (XI,*) ® 4>® (X 2 ,t), (7) 

where the coefficients C/ m , symmetrize the wave func- 
tion, and ® is the tensor product. The modulus squared 
of the coefficients |C; m | 2 , gives the probability of the pho- 
tons being in the states labeled by I and m. Each tensor 



component is related to the two-photon spin state's en- 
ergy probability density. The basis states j, are 

solutions of the single-photon wave equations (Eq. @ in 
free space and Eq.([4]) in a linear medium), and include 
spin dependence. The equation of motion for the two- 
photon wave function is found by adding the Hamiltoni- 
ans for the individual photons 

ihdM 2) = hv ia f ] (Vi + V1L1) x *( 2 ) 

+ hv 2 a% ) (V 2 + V 2 L 2 ) x *C9, (8) 

(2) (2) 
where a\ = £3 ® 1, a 2 = 1 ® £3, the curl operators 

are understood to act on appropriate components of the 

tensor product, L 1(2) = L(x 1(2 )), and w 1(2 ) = t>(x 1(2 )). 

In free space the Ls drop out and the speed of light takes 

on its vacuum value c. We call equation (|5J) the Maxwell- 

Dirac equation for a two-photon state. The two-photon 

wave function also obeys the divergence conditions 



(V, 



<J (2) = 0. 
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Tracing over the tensor product of the two-photon wave 
function and its Hermitian conjugate, and integrating 
over all space gives the expectation value of the prod- 
uct of the two photons' energies 



Tr 



\jr( 2 )t\|/(2) 



d 3 x 1 d 3 x 2 = (EiE 2 



(10) 



If the state of the photons is not entangled, this 
equals (E1E2) = (Ei) (E 2 ). One can also de- 
fine a joint probability density p^ (xi, x 2 , t) = 

Tr [tfW (x 1 ,x 2 ,t) t *( 2 3 (xi,x 2 ,t)l /(E 1 E 2 ), for finding 

the energy of one photon at the space-time coordinate x\ , 
and the other at x 2 , (xj = (x^ t), j = 1, 2), and current 
density (xi,x 2 ,i), obeying a continuity equation. 

To demonstrate the Birula-Sipe single-photon theory 
PI B El is best suited to PWM, we first show there is a 
direct relation between the n-photon wave function and 
the n-photon detection amplitude of quantum optics [j| 

E3. EE EE US 
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,x n ;i) = (vac\ <g> E (+) (x^t) 



\fr(«) 



(11) 

whose modulus-squared is proportional to the probability 
for joint, n-event detection. If one neglects the magnetic 
field, assuming the detectors respond only to the electric 
field, the effective ri-photon wave function is just a tensor 
product of n electric- field vectors evaluated at potentially 
different spatial values, exactly the same form as Eq. lfTTj) . 
Then the n-photon wave function can be identified with 
the spatial mode of the electromagnetic field, showing 
the connection between modes and states, and appealing 
to the choice of the single-photon formalism based on 
energy localization. 

To further strengthen the case for the energy-density 
single-photon theory, we show the close connection of the 
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two-photon wave function and classical coherence theory. 
Considering two positive-helicity photons for simplicity, 
the two-photon wave function can be written as a sum of 
four terms 



f (2) (xi,x 2 ) = D (a*) / V2e(xi) + iB (xi) /y/M^l 



D (x 2 ) /V2e(x 2 ) + zB (x 2 ) / VM^) ■ (12) 



Each term has the same form as one of the four second- 
order coherence matrices of classical coherence theory Q 



A (11,0:2) = (F* (xi)®G(x 2 )), 



(13) 



which give a complete description of second-order partial 
coherence of an optical field (including spatial, temporal 
and polarization coherence). Here F, G € {D,B}, and 
the brackets () imply an ensemble average over all real- 
izations of the fields. Evolution of these matrices is de- 
scribed by a set of linear differential equations @ that we 
call the first-order Wolf equations, which are equivalent 
to Eqs.© and ©. This equivalence shows a deep con- 
nection between propagation of classical coherence quan- 
tities and multi-photon states. In addition, each compo- 
nent of the coherence matrices obeys the (second-order) 
Wolf equations [9||, a well-known set of classical second- 
order differential equations recently highlighted for their 
relation to the two-photon detection amplitude [l^]. In 
much the same way that the Klein-Gordon equation does 
not completely describe the evolution of electron states 
by neglecting spin, the same holds for the second-order 
Wolf equations, which do not specify the relations be- 
tween polarization components. In contrast, the new 
Maxwcll-Dirac equation, Eq.®, contains all such rela- 
tionships. In this sense our result shows the quantum ori- 
gin of the classical Wolf equations [HI , and illuminates 
their connection to the propagation behavior of multi- 
photon states. Choice of a different single-photon theory 
upon which to base the multi-photon theory would not 
lead to the above results. 

To illustrate the utility of the two-photon wave func- 
tion and its relation to classical coherence theory, we con- 
sider the propagation through a turbulent atmosphere of 
two quasi-monochromatic photons, initially entangled in 
their spatial degrees of freedom, as depicted in Fig. 1. 
We assume the photons are emitted from a source in op- 
posite directions occupying one of two orbital angular 
momentum (OAM) states, described by the Laguerre- 
Gauss wave functions ip p j (r, 9) = (r) exp (U9) /V2n, 
where r and 9 are cylindrical coordinates. Here Rf (r) 
is the radial wave function, exp (U9) / ^/2~n is the angular 
wave function, I is the OAM quantum number, and p is 
the radial quantum number. We assume that both pho- 
tons, labeled A and B, have the same polarization, radial 
quantum number p — 0, and consider orbital quantum 
numbers of equal magnitudes |Z|, separately. We take the 
two-photon basis as 



<PAB 



v,-l 



p,—l ' 



(14) 



where tpfi 1S evaluated at the coordinate of photon 
A{B). For concreteness, we treat the input pure state 

^in = YPab+'Pab]/ 'v 7 ^- The photons pass through inde- 
pendent, thin, dielectric, Gaussian phase-randomizing at- 
mospheres, modeled by a quadratic phase structure func- 
tion [19l . |20| . This determines the form of the medium 
function I/(x) in Eq.([I)]), which we solve in the paraxial 
approximation. The elements of the density matrix p, at 

the output of the turbulence are determined by integrat- 

_ 2 

ing the radial power distributions r (r) , for each 

photon multiplied by the circular-harmonic transform of 
the phase correlation function C^, which describes the 
effect of the atmosphere on the state 21]. In our model, 
the phase correlation function for each atmosphere is 



cf B) (rA9) 



exp 



D 



MB) 



(r,A0) 



(15) 



where D^ (B) (r,A9) = [2r sin (AO) /r A{B) ] 2 is the 
quadratic phase structure function [20| of the aberrations 
in atmosphere A(B), and Oi(B) is the transverse length 
scale of the corresponding turbulence. We find a closed 
form expression for the circular harmonic transform 



cf B) (r, to) = C J l ya} (r, AO) exp {-imA9) dA9 



r 2n n A(B) 



exp 



"2 (r/r A{ 



B)) 



2 (r/r A{ 



where m is or 21, and I m (x) is an TOth-order modified 
Bessel function. 

Being interested only in two OAM states for each pho- 
ton, ip Pt ±i, we may treat each photon as a qubit. Other 
photon states can be considered to be loss channels [22j 
and we normalize the post-selected density matrix. We 
examine the decay of entanglement by calculating the 
concurrence C(p) [23| . from the normalized density ma- 
trix, as a function of the ratio of the optical beam waist 
to the characteristic turbulence length scale w/tq [2l| . 
Each atmosphere is assumed to have the same coherence 
length ro- For a maximally entangled state, C = 1, and 
for a non-entangled state, C = 0. We assume the atmo- 
sphere is unmonitored, so any independent information 
about its fluctuations is lost, leading to loss of entangle- 
ment. We plot the concurrence in Fig. 2(a), for various 
initial OAM quantum numbers. 

This result shows that for a beam waist much smaller 
than the turbulence length, w << ro, the entanglement 




/ = +l,±2,±3 




/ = +l,+ 2, + 3 



p,—l> 



FIG. 1: Diagram of a thought experiment involving two pho- 
tons, initially entangled in their OAM states, I — ±1, ±2, ±3, 
traveling through independent, random phase atmospheres, 
labeled by A and B. 
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FIG. 2: (Color online) (a) Concurrence as a function of the 
ratio, w/ro, of beam waist to turbulence length scale for three 
different magnitudes of input OAM. (b) Fidelity of the output 
with respect to the input pure state. 

is more robust to the turbulent atmosphere. Physically 
this reflects the fact that the photons will experience few 
phase distortions across their wave fronts. These results 
also indicate that entangled states with larger OAM val- 
ues experience less disentanglement through a turbulent 
atmosphere. This appears to be due to the fact that scat- 
tering from one OAM state to another depends only on 
the change in OAM, HAl [2l[, which must be supplied by 
the atmosphere. The atmosphere can, on average, change 
the OAM of the light only by a particular amount set by 
the spatial fluctuations that characterize it. We also cal- 
culate the fidelity of the output two-photon state relative 
to the input state, which for a pure-state input is 

F(pin,Pout) = (^in\Pout |*m) • (17) 

This result, plotted in Fig. 2(b) for the input state given 
above, indicates that states with small OAM values, and 



thus small "rms" beam width [2fj| . have higher overall 
transmission than do states with large OAM values. We 
should stress that the overall transmission of the OAM 
states depends on the beam waist w. We conclude that 
entangled states with smaller waists and larger OAM 
quantum numbers will be more robust to turbulence. 

We have introduced the two-photon wave function 
based on energy localization. This two-photon wave 
function obeys the two-photon Maxwell-Dirac equation, 
which is equivalent to the equations of motion of the 
classical second-order coherence matrices @. The con- 
nections we have found between this wave function, 
QED wave-packet-mode detection amplitudes, and clas- 
sical coherence theory give credence to the choice to use 
the energy-localization wave function rather than others, 
such as the Landau-Peierls wave function which is 
a non-local "number density" wave function. The for- 
malism provides powerful tools to analyze the behavior 
of few-photon states, as shown by the example above, 
where we calculated the disentanglement of a spatially 
entangled two-photon state by using essentially classi- 
cal field equations. This theory is well suited to the 
study of realistic implementations of linear optical quan- 
tum computing I24J . measurement-induced nonlinearities 
with linear optics [2a ] , and continuous- variable entangle- 
ment through quantum state tomography (2fj| . [27j . The 
well-defined Lorentz transformation properties of this 
wave function Q make it ideal for the examination of 
relativistic quantum information with photons [28| . 
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